The spin gap of CaV 4 Og revisited 
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The large-plaquette scenario of the spin gap in CaV4C>9 is investigated on the basis of extensive exact 
diagonalizations. We confirm the existence of a large-plaquette phase in a wide range of parameters, 
and we show that the most recent neutron scattering data actually require an intra-plaquette second 
neighbor exchange integral much larger than the inter-plaquette one, thus justifying the perturbative 
calculation used in the interpretation of the neutron scattering experiments. 
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CaV^Og is the only known two-dimensional magnet 
presenting a gap in the spin excitations. The low tem- 
perature dependence of the magnetic susceptibility first 
measured by Taniguchi et al. [Q provided a value of 
A ks 107 K. Since then several scenarios based on the 
Heisenberg model on the depleted lattice have been pro- 
posed to explain both the presence and the value of A. 
It appears that the presence of frustration induced by a 
second neighbor antiferromagnetic interaction is neces- 
sary and it has become clear that two types of first and 
second neighbors exchange integrals are required to fit 
experiments correctly. The Hamiltonian then reads, 



H = J\ ^2 Si.Sj + J[ ^2 Si-Sj 
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where the sums run respectively over bonds within pla- 
quettes (p) and dimer between plaquettes (d) for first 
and second nearest neighbor. 

The first explanation proposed to explain the pres- 
ence of a gap was based on a J\ — J 2 model on the 
depleted lattice with J2/J1 = 1/2 leading to weakly 
coupled small plaquettes. It must be abandoned for at 
least two reasons: on one hand, recent estimates of the 
exchange integrals in other members of the vanadates 
family suggest that the ratio J2/J1 is much larger than 
1/2. On the other hand, the most recent data obtained 
with inelastic neutron scattering seem to be inconsistent 
with the picture of weakly coupled plaquettes built on 
the smallest squares of the depleted lattice. In particu- 
lar, the dispersion curve of the lowest triplet state mea- 
sured by Kodama et al. Q turns out to be minimum at 
(0, 0) momentum. The authors then performed fits of the 
curve assuming that one of the integral exchange is large 
enough to justify a second order perturbation theory in 
the three other parameters and concluded that the best 
set of exchange integral is the one corresponding to J\, 
J[, J 2 small compared with J 2 , namely J\ = J[ 67 K, 



J 2 ~ 171 K and J 2 ~ 15 K. These parameters are also 
consistent with the temperature dependence of the sus- 
ceptibility if one assumes that part of the sample is non- 
magnetic H . In this picture the system consists of a set 
of weakly coupled J2-plaquettes. Nevertheless, regarding 
the mechanism responsible for the values of J 2 and J 2 
(super-exchange between two second neighbor Vanadium 
atoms via an oxygen atom) , the fact that these integrals 
differ by one order of magnitude is a priori surprising, al- 
though it is not necessarily contradictory since exchange 
integrals are very difficult to evaluate directly (from ab 
initio method in quantum chemistry for example). 
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FIG. 1. Symmetry of the first excited state for J[ — Ji. 
Symmetry sectors are labeled by the momentum k = (k x , k y ) 
and the irreducible representation of the axial point group. 
(0,0) A (S = 0) : (.), (0,0) B (S = 0) : (■), 
(0,0) B (S = 1) : (□), (0,0) x-iy (5 = 1): (V), 
(0,tt) A (S = 1) : (+), (tt.tt) A (S = 1) : (o), 
(ir,n) x + iy (S — 1) : (a). Note that the ground state 
is in the • sector. 



To assess the validity of the hypotheses underlying Ko- 
dama et aZ.'s analysis beyond perturbation theory, we 
have performed exact diagonalizations on a 16 site clus- 
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ter with periodic boundary conditions without making 
a priori assumptions on the value of the exchange in- 
tegrals. First of all, we have determined the symmetry 
of the ground state and of the first triplet excitation in 
the ( J 2/ Ji, J'il Ji) plane for different values of J[/Ji (see 
Fig. [I] for J[ = Ji). Surprisingly enough, it turns out 
that only rather small values of J 2 /Ji are compatible 
with the location of the minimum of the triplet band at 
k = (0,0). To check the robustness of this result with 
respect to variations of the ratio J[/J±, we have done 
a similar calculation as a function of (</(/ J±, J^/Jt) for 
J2/J1 = 2 (see Fig. ||). It shows that the restrictions 
are qualitatively the same, and that J 2 /Ji must always 
be much smaller than J2/J1 in order to reproduce the 
neutron scattering results. 
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FIG. 2. Same as Fig. 1 for J2/J1 = 2. 
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FIG. 3. "Phase diagram" in the (J2/J1, J2/J1) plane. The 
symbols denote points where (Sj.Sj) < —0.4 for 3\ (»),J[ 
(o),J 2 (■) and J' 2 (□) bonds. 



In conclusion, the results obtained by a systematic 
study of the energy spectrum and of the spin-spin cor- 
relation functions of a 16 site cluster in a large range of 
J'xl J\,Jil J\ and J' 2 j ' J\ ratios show that another gapped 
phase made of weakly coupled J 2 — plaquettes indeed ex- 
ists when J2/J1 is large enough. Besides, the experimen- 
tal fact that the minimum of the triplet band is located at 
k = (0, 0) implies that J 2 is much smaller than J2. This 
suggests that the hypotheses of Kodama et al. are indeed 
satisfied. Calculations on a 32 site cluster to confirm the 
present results are in progress. 
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Next, we have checked whether the set of parameters 
used in Ref. [0] and || was consistent with the picture 
of weakly coupled J 2 — plaquettes. To answer this ques- 
tion we have calculated the spin-spin correlation func- 
tion (Sj.Sj) with J{ = Ji in the presumed range of val- 
ues of J2/J1 and J 2 /Ji- Fig. || summarizes the results 
and clearly delimits three different regions in which the 
system is organized into small J\ -plaquettes, large J 2 - 
plaquettes and J[ — J 2 chains. We also checked that for 
J2/J1 > 2 the correlation functions on all bonds except 
J 2 never exceed 0.1 in absolute value, which shows that 
the coupling between J2-plaquettes is indeed very small 
in that range. 
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